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Boson-fermion mixtures realized in recent years with 
ultra-cold atomic gases have been known as a new 
many-body system exhibiting exciting phenomena and 
have attracted attention from both experimentalists and 
theorists 1 -. 4 The mixture systems are constructed not 
only in three-dimensional geometry but also in one and 
two dimensions. Several interesting features have been 
pointed out theoretically on the one-dimensional (ID) 
boson-fermion mixtures: Demixing instability 5,6 density 
wave, 7,8 Mott transition on a lattice, 9,10 and exact solu- 
tions for ground state properties. 11, 12 

In this short note we shed light on the mixing-demixing 
transition of ID boson-fermion mixtures. We could ex- 
pect the fermions and the bosons are demixed with 
strong-enough repulsive interactions between them. Ac- 
cording to a mean-field calculation, 5 the transition to the 
demixing phase occurs in a system with a low fermion 
density. On the other hand, a Luttinger liquid analy- 
sis revealed the possibility of the mixing-demixing tran- 
sition when the number densities of the fermions and 
bosons were largely different. 6 It was also shown in Bethe 
ansatz solution that the mixture was always in the mix- 
ing phase as far as the fermion-boson interaction had the 
same magnitude as the boson-boson interaction. 11 

In our previous work we numerically investigated the 
mixing-demixing transition of the boson-fermion mixture 
on a ID lattice at an incommensurate filling, 13 and found 
the transition always occurred when the fermion density 
was equal to or exceeded the boson density. 

Our aim here is to report the similar calculation in the 
other case, namely the case that the fermion density is 
below the boson density, extending the phase diagram 
obtained in Ref. 13 toward the negative side of dp — 
(Nf — N b )/N, where Nt, N b , and N are the numbers of 
the fermions, the bosons and the lattice sites respectively. 

We consider the mixture of Nf spinless fermions and 
N b bosons on an iV-site lattice by employing a Bose- 
Fermi Hubbard Hamiltonian: 

H = ^ [*o( a a,i ffl a,i+l + ' l - c -)+MQt n a 1 
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Fig. 1. The snapshot of the boson configuration (upper) and the 
fermion configuration (lower) with 5p = —4/30 and Ujt = 0.4 
(a) and 4 (b). X is the real space axis and r is the imaginary 
time axis. shows the number of the bosons (fermions) on 

each site. 
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where a a i and a a ^ are respectively creation and annihi- 
lation operators for fermions (a = /) or bosons {a = b) 
on the i-th site, and n a ,i = i(i a ,i- Hopping energy and 
chemical potential are denoted respectively by t a and p a . 
Ubb and Ufb are on-site boson-boson and fermion-boson 
repulsive interactions respectively. We set tf = t b = 1 as 
an energy unit. 

Reference 11 proved that the system always stays in 
the mixing phase as far as the fermion-boson interac- 
tion has the same magnitude as the boson-boson in- 
teraction, and also we reported in our previous work 13 
that the transition to the demixing phase occurred when 
the on-site boson-boson interaction is sufficiently smaller 
than the on-site fermion-boson interaction. In the present 
work, we thus set Ubb = to focus on how the transition 
point shifts with the change of dp. 

To observe the mixing-demixing transition, we mea- 
sured a correlation function defined by, 13 

C = {{n b .i-i+nb,i + i)n b .i) - ((n/.i-i+n/.j+n/^+iJn&.t). 

(2) 

When the two components are spatially mixed, the cor- 
relation function C tends to be smaller, or negatively 
larger, while it has a large positive value when they are 
demixed. Although we calculated other physical quanti- 
ties, C exhibited the change at the transition point most 
clearly and we decided to use C to locate the transition 
point. 

We performed Monte Carlo simulations with the world 
line algorithm on a ID periodic lattice. Temperature was 
fixed to T = 0.08, the Trotter decomposition number L 
to 100 and the number of sites ./V to 30. In the simula- 
tions, we changed the number of the fermions and the 
bosons, Nf and N b respectively, while fixing the total 
number Nf + N b = 14. As stated in the above we are 
interested in the case Nf < N b . 

At first, we show the snapshot of typical configurations 
for two components in Fig. 1, where the upper figures 
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area. Although we did not identify clearly the position 
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Fig. 2. The correlation function C as a function of Ufb in the 
cases of N f < N b (Sp < 0) (a) and N f > N b (Sp > 0) (b). The 
result in (b) was obtained in our previous paper 3 and used here 
for comparison. 



show the boson configurations and the lower show the 
fcrmion configurations when Sp = —4/30 and Ufb = 0.4 
in (a) and 4 in (b). With rather weak fermion-boson in- 
teractions, the system stays in the mixing phase as shown 
in (a) , while the strong interactions drive the system into 
the demixing phase as in (b). 

The transition can be seen more clearly in viewing 
the correlation function C, which is shown in Fig. 2 as 
a function of the fermion-boson interaction Ufb- Figure 
2(a) presents the correlation function behavior at Sp — 
-2/30,-4/30,-6/30, and -8/30. For comparison, we 
also show the previous result at Sp — 0, 2/30, 4/30, 6/30, 
and 8/30. As seen in the figure, there is a characteristic 
uprise in every curve of (a) and (b). This feature can 
be regarded as the occurrence of the transition from the 
mixing phase to the demixing at not only Sp > but 
also Sp < 0. It is interesting to note that the transition 
point seems to shift to the smaller Ufb side as the num- 
ber density of the fermions decreased relatively against 
that of the bosons. (Remember that the total number 
of the fermions and bosons was fixed in the simulation.) 
In other words, the shift of the transition point is not 
symmetric between the positive and the negative sides 
of Sp, but it is rather monotonic as Sp changes from the 
positive to the negative value. 

Our final aim of this short note is to complete the 
phase diagram in terms of the mixing and demixing 
phases. The phases could be determined by checking the 
snapshots and the behavior of the correlation function C, 
although the somewhat broad transition curves of C did 
not allow us to pinpoint the transition point. We there- 
fore classified each parameter point into three states: 
Mixing, demixing, and undetermined. The last one corre- 
sponds to the uprise part of the correlation function, i.e. 
the transition area (not "point"). The broad transition 
curves are caused by the finite size effect of the system, 
which we checked by changing the system size. 

Fig. 3 shows the phase diagram in the Sp — U fb plane, 
combined with the one in the area Nf > Nb obtained in 
the previous work. 13 The triangle symbols and square 
symbols denote the mixing states and the demixing 
states, respectively. The asterisks represent the transition 
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Fig. 3. Phase diagram in Nf — N b vs Ufb- The triangle sym- 
bols and the square symbols present the mixing states and the 
demixng states, respectively. The asterisks are marked for unde- 
termined phase points. The points of Sp > region were obtained 
in our previous paper. 13 The dotted line is just the guides for 
eyes. 



of the transition point, the figure strongly suggests that 
the transition point would shift to the smaller U fb side 
as Sp changes from the positive to the negative value. 
This behavior can be understood in the following way. 
Making a cluster in the demixing state is unfavorable in 
terms of the kinetic energy of the fermions, because the 
fermions cannot move around in the dense cluster. There- 
fore the demixing state requires sufficiently large Ufb to 
compensate the loss in the kinetic energy. However, as 
Sp changes from the positive to the negative value, i.e. 
as the number of the fermions decreases, the loss in the 
fermion kinetic energy becomes smaller, which makes it 
easier for the system to undergo the transition to the 
demixing state. This explains the shift of the transition 
point. 

In conclusion, we numerically studied the mixing- 
demixing transition of ID boson-fermion mixtures at low 
fermion densities with the total particle density fixed, 
and obtained the phase diagram in the Sp — Ufb plane. 
The phase diagram suggests that the qualitative feature 
of the system is not determined simply by the magni- 
tude of the difference in the number density between the 
fermions and the bosons. 

The authors would like to thank Prof. Jo and Prof. 
Oguchi for their continuous support. 
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